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Abstract Binary Systems are the most studied sources 
of gravitational waves. The mechanisms of emission and 
the behavior of the orbital parameters are well known 
and can be written in analytic form in several cases. 
Besides, the strongest indication of the existence of 
gravitational waves has arisen from the observation of 
binary systems. On the other hand, when the detec¬ 
tion of gravitational radiation becomes a reality, one 
of the observed pattern of the signals will be proba¬ 
bly of stochastic background nature, which are charac¬ 
terized by a superposition of signals emitted by many 
sources around the universe. Our aim here is to develop 
an alternative method of calculating such backgrounds 
emitted by cosmological compact binary systems dur¬ 
ing their periodic or quasiperiodic phases. We use an 
analogy with a problem of Statistical Mechanics in or¬ 
der to perform this sum as well as taking into account 
the temporal variation of the orbital parameters of the 
systems. Such a kind of background is of particular im¬ 
portance since it could well form an important fore¬ 
ground for the planned gravitational wave interferome¬ 
ters DECI-Hertz Interferometer Gravitational wave Ob¬ 
servatory (DECIGO), Big Bang Observer (BBO), Laser 
Interferometer Space Antenna (LISA) or Evolved LISA 
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1 Introduction 

Of all the theoretically possible sources of gravitational 
waves, the compact binary systems, that is, those sys¬ 
tems whose components are black holes or neutron stars, 
are the most interesting, both from the theoretical and 
observational viewpoints. Goncerning the theoretical as¬ 
pects, it is possible to derive equations which hold even 
in the advanced stages of spiralling, while from the ob¬ 
servational viewpoint, binary systems are among the 
most probable sources to have their gravitational radia¬ 
tion detected in the near future by the planned gravita¬ 
tional wave interferometers DECI-Hertz Interferometer 
Gravitational wave Observatory (DECIGO), Big Bang 
Observer (BBO), Laser Interferometer Space Antenna 
(LISA) or Evolved LISA (eLISA), Advanced Laser In¬ 
terferometer Gravitational-Wave Observatory (ALIGO) 
(and VIRGO), and ET, as we will see later on. 

On the other hand, concerning the study of the grav¬ 
itational waves themselves, of special interest for us are 
the stochastic backgrounds, which are characterized by 
spectra spanning a wide range of frequencies, where we 
cannot distinguish individual sources, instead, we can 
distinguish only a “smooth” shape, characteristic for 
each type of physical process which created them. 

Here we are interested in the background of grav¬ 
itational waves generated by compact binary systems 
during their periodic or quasiperiodic phases in a cos¬ 
mological scenario. The background considered here is 
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similar to the confusion noise generated by Galactic bi¬ 
naries in the LISA or eLISA frequency band[T]. 

Before proceeding, it is worth mentioning the well- 
known fact that other authors studied such a subject. 
For example, in a recent review paper, Regimbau[2] 
considered the astrophysical gravitational wave stochas¬ 
tic background, where the background by the compact 
binaries is also considered. Recently, Zhu et alia consid¬ 
ered the stochastic gravitational wave background from 
coalescing binary black holes and Rosado [3] the grav¬ 
itational wave background from compact binary sys¬ 
tems in general. It is also worth mentioning some earlier 
works on this subject, including Farmer and Phinney[a, 
Regimbau and de Freitas Pacheco 0 and Schneider et 
aim. Some ingredients considered by many of these au¬ 
thors include, for example, population synthesis and 
cosmic star formation histories. 

Given that many authors studied the background 
generated by compact binaries, one could ask whether 
it is still worth studying such a subject. Since there is 
other possible approaches, such as the one considered 
in the present paper, there is still space for keeping 
studying such a subject. 

Then, we shall present here an alternative method 
to calculate such a background, which is organized in 
steps. First, in Section[21 we present some properties of 
the gravitational radiation emitted by binary systems 
such as frequencies, intensities and their relation with 
the orbital parameters; in Section [31 we describe the 
population characteristics of the systems in the form 
of distribution functions involving orbital parameters 
such as orbital distances and masses of the components 
of the systems and the star formation rate. Once we 
have individual and population characteristics, we are 
ready to calculate the stochastic background, and in 
Sect, m a method for this calculation is shown, where 
we adapted a method used in the solution of a simple 
problem of Statistical Mechanics in order to obtain the 
desired background spectrum; in Section [5l we present 
and discuss some results obtained by using the method 
mentioned above, besides showing and justifying the 
numerical parameters we used in the calculation; in Sec¬ 
tion El we consider the detectability of the spectra by 
the cross-correlation of pairs of LIGO and ET, and in 
Section [71 we present some conclusions and perspec¬ 
tives. Finally, in order to give the reader more details, 
there is one Appendix at the end of the paper. 


2 Binary Systems and Gravitational Waves 

In the present paper, we are concerned with systems in 
circular orbits, where the only harmonic present in the 
emitted radiation is the first one, whose frequency is 


twice the orbital frequency; for non null eccentricities, 
the higher harmonics will be in general important, in 
such a way that a fraction of the total power will be 
emitted at higher harmonics [SllQlfTO] . 

For circular orbits, the expression for the total power 
emitted in gravitational waves by binary systems is 
given by mi: 

dE 32GI?® / 17111712 

dt 5c® + *^2 

where mi and m 2 are the masses of the components of 
the system, r is the orbital distance and 17 is the orbital 
angular frequency which can be written as a function 
of the orbital distance via Kepler’s third law: 

17 ^ = G{mi+m2) ^ ^ 2 ) 

Besides, it is possible to put 17 as a function of time; 
with the aid of o and after some algebraic manipula¬ 
tion, we obtain 

( 3 ) 

dt 5 c® 

whose integration yields the orbital frequency as a func¬ 
tion of time, namely, 

= (4) 

In ([dj, 17o is the initial frequency, to is the initial 
time and K is a constant given by 

5c® 

This expression, even though valid only for circular 
orbits, will be useful because the method developed in 
the present paper for the calculation of the stochastic 
background was first applied to circular orbits. Besides, 
it gives us a good physical insight on how the evolution 
of the systems takes place. 



3 Population Characteristics 

Binary systems that contribute to the background may 
present orbital parameters such as semimajor axis and 
masses of the components that are quite different from 
system to system, which generally fill a continuous range 
of values. In general, it is possible to define mathemat¬ 
ical functions which describe how the several values of 
a given parameter are distributed among the elements 
of the population. Such distribution functions are vital 
to the calculation of the background and may be found 
in the literature under several forms. 
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3.1 Mass Distribution Function and the Star 
Formation Rate 


index n = 1 , normalized to the abundance of rich galaxy 
clusters at present day (erg = 0.9) - found a star for¬ 
mation rate covering a range of redshifts from 2 ; = 0 to 

We are considering here the so called Salpeter distribution fl^ p “dark ages” {z ~ 20). Such a function is given by 

P^a(z-z^) 


which is written as 


( 6 ) 


P*{Z) = Pr, 


(9) 


where x = 1.7 and A = 0.17, with obeying the 

following normalization condition 

r™/ 

/ m(j){m)dm = 1. (7) 

J rrii 

Here, m/ and rrii are the maximum and minimum val¬ 
ues of the masses that stars can have, respectively. 

We consider that stars can have minimum masses 
of O.IMq and maximum masses of 125M0 where 
this interval is split in subdivisions: stars with masses 
between O.IMq and 8Mq give rise to white dwarfs at 
the final stage of their evolution; stars with masses in 
the interval between 8Mq and 25Mq originate neutron 
stars and stars more massive than 25Mq end their ex¬ 
istences as black holes. 

It is worth mentioning that Salpeter distribution is 
valid only for the masses of the progenitor stars, and 
when such progenitors become black holes or neutron 
stars, their masses change, so we must estimate the val¬ 
ues of the masses of such remnants. For the mass of the 
neutron stars, we use a constant value of 1.4 Mq, ac¬ 
cording to the Chandrasekhar limit and the theories of 
formation of supernovae II and Ib/c|13j. For the rem¬ 
nant black hole mass distribution, we use found in 
a recent study by Ozel et alia, namely 


Gbh(mbh) = 0.332exp [-0.347(TObh - 7.8)^] , 


( 8 ) 


_ Q, -P 

where a = 3/5, j3 = 14/15, Zm = 5.4, and with pm = 
0.15MQyr“^Mpc“^ fixing the normalization. Further, 
it is worth mentioning that the authors used hydrody¬ 
namic simulations, taking into account radiative heat¬ 
ing and cooling of gas, supernova feedbacks and galactic 
winds. 

Even though ([5]) gives information about the star 
formation up to 2 = 20, we must point out that the 
observations can give information from at most 2^5. 
Therefore, uncertainties are inherent in any modeling 
above a given redshift. Also, the star formation rate 
adopted in this work refers to the first population II 
stars formed; therefore it is consistent with the Salpeter 
IMF adopted. 

On the other hand, in order to see how the spectra 
depend on this important ingredient, we also performed 
calculations using other SFRDs. In particular, we con¬ 
sidered four additional star formation rates [IH]. Three 
of them are given by 

p*[z)i = l.Q^C^hQ^F{z)G^{z)MQyv~^Mpc~^ (10) 

with i = 1,2,3 denoting each case, = /i/0.65, h = 
0.7, F{z) = [Cm(1 + z)^ + Ci = 0.3, 

C 2 = 0.15, C 3 = 0.2, Ga = 0.7, Cm = 0.3 and the 
functions Gi{z) given by 

„3.4z 


^ 3.82 


45 


c3.4z 


where the masses are given in solar masses. Concerning 
the range of masses of black holes, we consider the range 
of 5Mo < mbh < 20M^ fT4lfT5| . However, we must em¬ 
phasize that there are uncertainties in the determina¬ 
tion of exact values of this interval of masses, so we 
chose values that seem suitable, bearing in mind the 
cited references. An interesting study concerning the 
black hole formation parameters can also be found in 
Postnov and Yungelson[lfij. Later on, we comment on 
how our results would be affected if a wider range of 
values were assumed. 

However, the description of the star population is 
not complete without the star formation rate density 
(SFRD), and we use the function by Springel and Hernquist 
where the authors - considering a ACDM cosmology 
in a structure formation scenario with the following 
cosmological parameters: Cm = 0.3, Ga = 0.7, Hub¬ 
ble constant Hq = 100/ikms“^Mpc“^ with h = 0.7, 

Cb = 0.04 and a scale invariant power spectrum with 


a3.4z 


Gi{z) = 

G 2 {z) = 

= g2.93z 15 

The fourth one has the following form: 


22 


,^3.052-0.4 


p^{z) = h 


0.017-b0.13z 


MQyr ^Mpe 


..-3 


( 11 ) 

( 12 ) 

(13) 

(14) 


A-b (z/3.3)5-3 

It is worth mentioning that we will not consider here 
a complete description of these additional SFRD’s. We 
refer the reader to the paper by Zhu et al[TH] for a 
detailed discussion concerning this issue. 

Notice that the three first functions describe forma¬ 
tion rates for redshifts of up to z ss 4 and the fourth 
121^ is valid for redshifts of up to z « 6. Besides, we use 
the following notation[T8]: for the rates given by (dQD 
we use the labels SFRl, SFR2 and SFR3 and for (fT4l) 
we use HB06; the formation rate by Springel and Hern¬ 
quist is written as SH03. Figure [T] shows the SFRDs 
adopted in the present paper. 
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redshift z 

Fig. 1 The SFRD’s adopted in the present paper. See the 
text for details 


Table 1 Parameters of the distribution functions of the or¬ 
bital distances 


system 

C 

r(R0) 

cr(Ro) 

DNS 

0.070 

0.6 

0.2 

BHNS 

0.015 

5.5 

1.5 

BBH 

0.070 

11.0 

2.5 


For black hole-neutron star (BHNS) and binary of 
black holes (BBH) systems, we consider that the pop¬ 
ulation of compact binaries is composed by DNS sys¬ 
tems (61%), BBH systems (30%) and BHNS systems 
(9%) [in]. Therefore, we may use these proportions to 
estimate the mass fraction for BHNS and BBH systems. 


3.2 The Binary Formation Rates 


3.3 The Semimajor Axis Distributions 


Before proceeding, it is necessary to know the rate of 
formation of compact binary objects; that is, given a 
population of progenitor stars, we need to know the 
fraction of stars which generate double compact bina¬ 
ries. As we may suppose, such fractions may assume 
different values because, from the formation of a popu¬ 
lation of stars to the final configuration (which is char¬ 
acterized by a population containing main sequence 
stars, neutron stars, black holes and binaries formed 
from these objects), several processes may affect the fi¬ 
nal number of compact binaries present. For example, 
if a component of a binary system gives rise to a su¬ 
pernova, the explosion may break the system. Further, 
note that not all progenitor stars form binary systems, 
and not all stars will generate black holes or neutron 
stars. 

For the formation rate of double neutron star (DNS) 
systems, we consider the mass fraction Adns that is con¬ 
verted into DNS systems [B], which is given by 

Adns — Pnsfp^ns (1^) 

where /?ns is the fraction of binary systems that survive 
to the second supernova event, fp gives us the fraction 
of massive binaries (that is, binary systems where both 
components could generate a supernova event) formed 
inside the whole population of stars, and ^ns is the mass 
fraction of neutron star progenitors that, in our case 
and using ([B]), is given by 

1-25 

^ns = / 4 >{m)dm. (16) 

Js, 

Numerically, we have /3ns = 0.024, fp = 0.136 and <?ns = 
5.97 X lO-^Mg^ 

Finally, the binary formation rate for DNS systems 
is given by 

^dns('2') = AdnsP*(^)- (1^) 


For the distribution function of the semimajor axis, we 
followed Belczynski et al|19). who made a very compre¬ 
hensive study of the compact binary systems. In partic¬ 
ular, we used the results present in Fig. 4 of that paper. 
To deal more easily with the distribution functions ob¬ 
tained by Belczynski et al, we ht them with gaussian 
functions. A close comparison between our functions 
(see Fig. [2|) and Fig. 4 of the cited paper shows that 
the fits work very well. 

The semimajor axis distributions are then written 
as follows: 


/(r) = Cexp 


'-(r - r)^' 

2 ^^ 


(18) 


where r is the semimajor axis, which is given in solar 
radius units, and the parameters C, f and tr are given 
in Tabled] 

It is worth noting that the above distribution func¬ 
tions obey the following normalization condition|19|: 


[fdnsir) + /bhns(c) + fhhhir)]dr = 1 


(19) 


In addition, it is useful to have information about 
the orbital periods or frequencies of the systems. A con¬ 
venient way to achieve this is to plot the distribution of 
the initial frequencies. As separations and frequencies 
are closely related by Kepler’s third law, it is straight¬ 
forward to obtain the distribution of the frequencies 
from (flSl) (see the next Subsection). 


3.4 Time Evolution of f{r) 

A remarkable feature of binary systems is the time evo¬ 
lution of the orbital distance due to the emission of 
gravitational radiation. For the systems considered in 
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Fig. 2 Distribution functions of orbital separations of com- Fig. 3 Distribution g{u) for NSNS systems 
pact binaries 


this paper, such variations are very important, espe¬ 
cially at the stages near the coalescence, so the dis¬ 
tribution functions shown above must be modified in 
order to take into account the time as an independent 
variable. First, let us consider the distribution given by 
®, which will be modified with the aid of (|3]), valid 
for circular orbits. 

Before including the time variable in (ITR)) . it is con¬ 
venient to put this function in terms of the frequency, 
by means of ©■ So, changing the variables via 

f{r)dr = g{f2)df2 (20) 

we get after some algebra 


ff(^) 


T 


X [G(mi -I- 1712 )]^^^ Q ®'^^exp 


— {r — r) 


21 


2a2 


( 21 ) 


Now we use and perform a change of variables 
in (EU, writing 

g(l7o)dl7o = H{Q)dn (22) 


where l7o is the initial frequency, which was associated 
with the variable 1? in (ED). So, performing the deriva- 
tive of (El) we have 


dfl \ ^ ) 


(23) 


which, substituting in (1^ . gives after some algebra 


H{n) = ^ [G{mi + 1712 )]^^^ 


xl7 ^^/^l7Qexp 


-{r - r)2 


where 

= 


G-®/3 + -Kit - to) 


2a2 

-3/8 


(24) 


(25) 


and K is given by m- We should do a further coordi¬ 
nate transformation, writing H[f2) as a function of the 
emitted frequency v instead of the angular frequency 
1?. Such a transformation, given by z/ = l?/7r, is trivial, 
and all the equations used from now on will be writ¬ 
ten as a function of v. As an example. Fig. [3] shows the 
initial distribution g{v) for DNS systems. 

Note that r and f were maintained in the above 
equations just for simplifying their presentations; they 
are effectively substituted by the appropriate set of vari¬ 
able transformations. 

It is worth noting that Hils et al[2Q] followed a sim¬ 
ilar procedure; namely, they started with a given initial 
distribution function and evolved it as a function of 
time. 


3.5 Cosmological Parameters 


The cosmological model adopted here, which is the ACDM, 
has the volume element given by 

dz 

where 

1 


— = 47r ( — ) r^Fi^M, Ga, z), 


(26) 


F'(17m, z) = 


\/( 1 + 2 ;) 2 ( 1 -|- I^mz) — z(2 + z)Ga 

and the comoving distance given by 
c 


(27) 




\/|Ck| / Fif2M,GA,z')dz' 


(28) 


where the density parameters (which should not be con¬ 
fused with the frequency 17 in the previous sections) 
obey the following relations: 


17m = 17dm + 17b and 17 m + l7yi + l7k = 1 (29) 


The subscripts M, DM, B, A and k refer to matter, 
dark matter, baryonic matter, cosmological constant 
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and curvature, respectively. Furthermore, the function 
S is defined by 

{ sin(a;) closed universe 
X flat universe 
sinh(a;) open universe 

In the present paper we consider the case of a flat 
Universe (i.e. I7k = 0), where the cosmological param¬ 
eters were presented in the previous section. 

4 The Stochastic Background 

The stochastic background of gravitational waves gen¬ 
erated by compact binaries will be calculated by means 
of|21) 

^BG = ~ / (30) 

^obs J 

where /ibg represents the dimensionless amplitude of 
the spectrum, /isource is the amplitude of the emitted 
radiation of a single source, i^obs is the observed fre¬ 
quency, and dR is the differential rate of generation of 
gravitational radiation. The differential rate of produc¬ 
tion of gravitational waves is derived by using a sta¬ 
tistical approach. One of the aims of such a derivation 
was to take in to account the change in the gravita¬ 
tional wave frequency emitted by binary systems that 
are particularly important for the system emitting at 
the high band of the spectra. 

In the present case ft-source is given bv pilH] 

^source = 7.6 X 10“^^ 



where ^ is the reduced mass of the system, M is the 
total mass, and du is the luminosity distance. Besides, 
the differential rate dR in pOD for the case of BBH 
systems may be written as 

dR 

d7?bbh = —Gh\v{mi)Gh\i[m2)dzdmidm2 (32) 

dz 

Notice that we consider in the above equation the 
most general case, where the masses of the components 
of the systems are not constant, which in our case refers 
to the systems which are formed by two black holes 
obeying the distribution function given by (|H]). For the 
case where one or both components are neutron stars, 
one or both distributions for the masses in (1321) will be 
considered as a Dirac delta, since we are considering 
that all neutron stars have the same masses. Further, 
note that, in order to simplify the notation, we have 
written the masses of the black holes as mi and m 2 . 


In addition, the term dR/dz will be determined in a 
subsequent section. 

Still concerning (I30|) . it is worth mentioning that 
the starting point of its derivation comes from an en¬ 
ergy flux equation. This equation was first derived in a 
paper by de Araujo et al[53] and was used in their var¬ 
ious articles. In particular, in a subsequent paper pT| 
they gave a more detailed derivation of this equation, 
showing its robustness. Also, although apparently sim¬ 
ple, it contains the correct and necessary ingredients to 
calculate the background of gravitational waves for a 
given type of source. 

It is also important to call attention to a particu¬ 
lar issue considered in the cited paper [21]. They raised 
a discussion regarding the inclusion or not of an addi¬ 
tional (1 + .^) factor in the denominator of the integrand 
of (15(11) . which many authors consider to take into ac¬ 
count the time dilation due to the expansion of universe. 
These authors argued that when one uses the luminos¬ 
ity distance in the calculation of the gravitational wave 
amplitude ^source, this distance already contains all nec¬ 
essary (1 -(- z) factors. This is in fact a controversial 
issue. In practice, however, the inclusion of the addi¬ 
tional time dilation factor does not modify significantly 
the amplitude of the gravitational wave background. 
Therefore, in the present paper we included the (1 -I- z) 
factor in the denominator of the integrand of (1501) . 


4.1 The Statistical Problem 

Basically, the problem of determining the stochastic 
background comes down to the counting of the number 
of binary systems that, at a given instant of time, emit 
gravitational waves at a given frequency. Specifically, 
when one is dealing with high frequencies (about IHz or 
more), the systems evolve very fast, and the frequencies 
of the emitted radiation vary in a small timescale, mak¬ 
ing the calculation of the systems which have the same 
frequency at a given moment of time more difficult. 
However, the problem can be handled by adopting a 
particular way of modelling the behavior of the popula¬ 
tion of binary systems; let us idealize a one-dimensional 
“space” where the values of the coordinates are the val¬ 
ues of the emitted frequencies. 

In this context, the evolving binary systems can be 
represented as points moving in the direction of increas¬ 
ing values of the coordinate (emitted frequency). At this 
point, we can note that those points have a behavior 
similar to a gas inside a chamber, where the particles 
have a certain spatial distribution and move towards 
a given wall of the chamber at different velocities. So, 
using this analogy, we can adapt some techniques of 
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Fig. 4 The volume element dV adjacent to A 


statistical mechanics in order to develop a method to 
calculate the background. 

Let us consider a gas in thermal equilibrium con¬ 
tained in a rectangular chamber of cross section (rela¬ 
tively to the X axis) equal to A. An elementary ques¬ 
tion concerning such a configuration is to calculate how 
many particles (or molecules) reach a given surface (per¬ 
pendicular to the X axis) of area A in a time interval 
dt; that is, the aim is to calculate the flux F of particles 
through a given area. Discussions on this question may 
be found in the literature, and we follow a particular 
way of handling the problem|24): let then Adx = dV be 
a differential volume adjacent to the surface A, accord¬ 
ing to Fig. m 

If we consider that the particles are distributed in 
a homogeneous way inside the chamber, the fraction of 
particles inside the volume dV is given by 



Adx 


(33) 


where V is the total volume of the chamber, but using 
the velocity of a given particle at the x direction we 
have: 


dV = Avxdt (34) 

In this way, all the particles inside the volume dV 
will reach the surface A in a time interval dt. In order to 
calculate the total number of particles that hit the area 
A, we should integrate over the velocities Vx, using a 
given velocity distribution function, which we write as 
'r]{v). So, considering that r]{v)dv^ is the total number 
of particles at a given interval dv^, we get 

d/i = —— r]{v)d^v (35) 

where d/r is the number of particles in dV, which have 
velocities in the interval d^v. Now, we should integrate 
only over the positive values of Vx, because we are con¬ 
sidering only particles that are moving towards the sur¬ 
face A, so we get the number of particles per area A and 
time interval dt (the desired flux F): 

^ pOO pCO pCO 

F = — dvx dvy dv^Vxri{v) (36) 

^ Jo J—oo J— cx> 


This is sufficient for the solution of the statistical 
problem, but we need to make two modifications in the 
above derivation in order to adapt it to the calcula¬ 
tion of the stochastic background. First, note that the 
cross section of the chamber is constant, so dMl) may 
be rewritten as 


dn = 


Adx 

~AL 


dx 

T 


(37) 


where L is the total length of the chamber. Thus (1551) 
is rewritten as 


dfi 



(38) 


A second modification is to consider a more general 
case, where the spatial distribution is not constant. So, 
it is reasonable to write the expression above as 

dfi = ip{x)r]{v)d^vdx (39) 


where il){x)dx is the fraction of particles at the interval 
dx. If the spatial distribution obeys a function which 
depends only on x, which we call ip{x), we can write 
the function tplx) in the following way 


iIj{x) 


‘Pjx) 

J (fi(x)dx 


(40) 


for the case where ip(x) is not normalized. With these 
two considerations at hand, the problem is ready to be 
used in the calculation of dR. 


4.2 Calculation of dR 


Before proceeding, it is worth stressing that the deriva¬ 
tion of dR presented here refers to the source frame. 
Therefore, ly in the equations below is the emitted fre¬ 
quency, and the time derivatives are also taken with 
respect to the source frame. 

The calculation of dR by means of the statistical ap¬ 
proach shown above is done by performing two changes: 
instead of using the spatial distribution <p(x), let us 
use a frequency distribution (p(iy) and change the ve¬ 
locity distribution for a distribution in a new variable 
Ujy, which we shall define by 


diy 

dt 


(41) 


that is, the time derivative of the frequency (in anal¬ 
ogy to the definition of velocity, which establishes the 
time variation of the spatial coordinate). Once this is 
done, we need to find out the forms of these two new 
functions. 

So, the function ip{iy) may be defined as follows: 


ip{iy) 


J nhin{to)H{t, to, iy)dto 


(42) 
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where to represents the instant of the birth of the sys¬ 
tems, nbin(^o) is the formation rate of binary systems 
(where the subindex bin is to be substituted by dns, 
bhns and bbh in each case) and H{t,to,v) is given by 
(123 and represents the fraction of systems that begins 
to exist at the instant to and that have frequencies in a 
given interval dv. Note that here, for sake of notation, 
we explicit the dependence of (IMl) on t and to- Besides, 
we refer the reader to Appendix |X] to see the derivation 
of (HU)- Further, note that Zo and to are related to each 
other by to = + (1 + 

The distribution r]{vi,) will have a peculiar form. 
First, we need to note that from (ED we have 

(43) 

Then, we conclude that there will be just one value of 
Vi, for each value of which allows us to write r]{vi,) 
as a Dirac delta function 


r]{v^) = N5{vu - Vv,p) (44) 

where N is the total number of systems and Vv,p 
is the particular value of Vi, corresponding to each fre¬ 
quency V. 

Now we can rewrite as 


dfj. = (45) 

where we used (HOD with the variable x substituted by v. 
Now, noting that the denominator of the term between 
parenthesis is the total number of systems, using the 
function given by (l4^ and changing the differential dv 
by means of the chain rule, we get 

dii= i ^ N5{vi, - Vv,p)dv^ (46) 


By integrating over Vi^ and rearranging the result, we 
obtain 




(47) 


where R is number or systems by time interval dt. Re¬ 
calling that this rate R is per comoving volume, we may 
write explicitly 


dv dR 
^ dV 

and (l30ll assumes the form 


^2 _ 


^obs 


"source 


dR dV 
dVl + z' 


(48) 


(49) 


Further, using this amplitude we can derive a largely 
used quantity, the so called spectral amplitude, which 
is given by: 

Sn = (50) 

^obs 


It is worth noting that we use the classical (i.e., non- 
relativistic) approach in the above statistical problem. 
There are two reasons for this choice: first, we should 
note that even in the late stages of evolution, the com¬ 
ponents of the binary systems have low orbital velocities 
(i.e., nonrelativistic velocities), and second, in all stages 
of evolution, the time variation of the orbital frequency 
(ui/) is not fast enough, and therefore, it is not nec¬ 
essary to consider special relativistic corrections when 
one makes the analogy between v^, and the velocity of 
the particles. Concerning general relativistic effects, the 
well-known quadrupole formula is more than enough. 
This formula is (implicitly) used in the calculation of 
the time evolution of the orbital frequency (see, e.g., 

HI)- 

Besides, we refer the reader to the paper [55], where 
we show the (first) simple example of application of 
the method described here; namely, we calculate the 
stochastic background generated by double neutron star 
systems, considering an uniform orbital period distribu¬ 
tion. 


5 Results and Discussions 

One of the aims of the present paper is to investigate 
whether the background of gravitational waves by the 
cosmological compact binary systems could in princi¬ 
ple generate a “confusion noise” for the spatial and ter¬ 
restrial interferometers; in the same way the Galactic 
binary systems do for LISA[T]. If such confusion noise 
really exists, it must necessarily be taken into account 
since it could well dominate the sensitivity curves of the 
detectors. 


5.1 Setting the Maximum and Minimum Frequencies 


For DNS systems, we consider for the maximum fre¬ 
quency of the gravitational waves the emitted value of 
« QOOHzEJ; for BHNS systems, we use (IHTll in order 
to determine the maximum frequency 

G(mns -f m-bh)’' 


isco 


(51) 


where rrins = 1.4M0, Wbh = 5.OM0 and risco is the in¬ 
nermost stable circular orbit (ISCO) of the black hole [55] 
given by 

^ (52) 

For BBH systems, the maximum frequency is given 
by|29] 


ao'n'^ 


boT] + Co 


G 


ttM 


( 53 ) 
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Spectral Amplitude for NS-NS Systems 



logv 

Fig. 5 Spectral amplitude for DNS systems and the sensi¬ 
tivity curves of some interferometric detectors 


Spectral Amplitude for BH-BH Systems 



logv 


Fig. 7 Spectral amplitude for BBH systems and the sensi¬ 
tivity curves of some interferometric detectors 


Spectral Amplitude for BH-NS Systems 



logv 


Fig. 6 Spectral amplitude for BHNS systems and the sensi¬ 
tivity curves of some interferometric detectors 

where 17 = mbhiWbha/-^^ is the symmetric mass ratio, 
M is the total mass and the polynomial coefficients are 
ao = 2.974 x 10"^ bo = 4.481 x and cq = 9.556 x 
10 - 2 . 

The values of the minimum frequencies do not play 
a decisive role in the calculation; that is, they do not 
affect the shapes of the spectra. Therefore, we consid¬ 
ered the value of 10 “® — 10 -®Hz for the three fami¬ 
lies of compact systems. This choice is suggested by 
analyzing the behavior of (l 2 T]l (see also Fig. |3l). As 
can be seen, the number of systems emitting below 
10 “® — 10 -®Hz is negligle as compared with those emit¬ 
ting around 10 “^ — 10 -®Hz. 

5.2 The spectra 

Figures [SJIH and [7] show the spectra for DNS, BHNS 
and BBH systems, respectively, which are compared 
to the sensitivity curves of LISA, eLISAjSO], BBO|3T], 
DECIGO[32], ET[33] and ALIGO[33]. The sensitivity 
curve for LISA may be found at 
WWW.srl.caltech.edu/~shane/sensitivity/ 


When observing the curves for the DNS and BBH 
systems, we note that LISA, eLISA, BBO, DECIGO, 
ALIGO and ET would not be affected by such a back¬ 
ground. BHNS binaries would be marginally important 
for ET at the frequencies of highest sensitivity. 

When observing the shape of the spectra, we note 
that there is a cutoff frequency, which in the present 
case is equal to the maximum frequency of the gravi¬ 
tational waves allowed for each type of system. Beyond 
this maximum frequency, the systems reach the chirp 
phase that is characterized by a particular pattern of 
emission, but we are not considering such a phase here. 
In fact, the chirp and coalescing phases will be consid¬ 
ered in forthcoming papers. 

It is also interesting to see the influence of the ini¬ 
tial distribution function of frequencies on the spectra. 
As the distributions for orbital separations and frequen¬ 
cies are closely related, we can achieve this objective by 
handling the distribution of the initial separations. So, 
as an example, we calculated the spectra for DNS sys¬ 
tems, considering the hypothetical cases where for the 
standard deviation a of the initial distribution given 
by (fT^ gives the values 2a and a 12. As expected, the 
semimajor axis distribution has little influence on the 
spectra for frequencies much shorter than the smallest 
typical birth period. 

Our calculations also show that the resulting spec¬ 
tra do not depend significantly on the SFRDs adopted. 
There are two main reasons for this result. The first 
one has to do with the fact that the sources located at 
2 > 5 do not significantly affect the resulting spectra. 
The second reason is related to the fact that the SFRDs 
for z < 5 are similar. 

Besides, through m, we observe the influence of 
the masses of the systems. For comparison, let us con¬ 
sider a BBH system where both components have the 
maximum value of mass, but considering two cases: in 
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the first case considering that the maximum mass is of 
^ 6 OM 0 [ 33 ], and in the second case, using the more 
realistic value of 20Mq as an upper limit (remembering 
that this is the value considered in the present paper). 
The amplitude given by (EID in the first case is 6 
greater than the amplitude given by the second case. 
As a result, the amplitude of the background would 
also be higher. 

Notice that the shapes of the spectra are different. 
The reason for this is the following. For a given semi¬ 
major axis, DNS systems emit gravitational waves at 
the same frequency, since we consider that all neutron 
stars have the same masses. For BBH (BHNS) systems, 
however, the emitted frequencies depend on the black 
holes masses that, contrary to the neutron stars, as¬ 
sume different values since they follow a distribution 
function. 

It is worth paying attention to some similar studies 
found in the literature, since they could well consider 
different approaches. For example, in Schneider at al[7] 
the spectra were calculated for binary systems formed 
by neutron stars, black holes and white dwarfs at the 
early stages of spiralling. These authors used a differ¬ 
ent method, based on binary population synthesis pro¬ 
grams. The differences in the results have to do with 
the fact that in their calculations, the number of sys¬ 
tems that contributes to the lowest frequencies is much 
greater than that of those that contribute to the high¬ 
est frequencies. Furthermore, these authors considered 
lower values for the minimum frequencies than those we 
are using. 

In Regimbau[35], the stochastic background gener¬ 
ated by DNS systems in the frequency band of LISA is 
considered, where the author used a particular proba¬ 
bility distribution[36) in order to generate a population 
of DNS systems. Here, the author considered the time 
evolution of the orbital parameters in a different way 
we do; actually, she considered the difference between 
the redshift of the birth of the systems Zh and the red- 
shift of the emission Ze of gravitational waves, where z^ 
and Ze are related to each other by the frequency of the 
emitted waves. 


6 Cross Correlation of Detectors 

The predicted spectra of the present paper cannot be 
detected by single interferometric detectors such as the 
ones cited above, but a putative detection of these spec¬ 
tra can occur by means of a suitable correlation of 
two or more detectors [37l[38l|39] . More specifically, we 
will consider the cross correlation for pairs of Initial 
LIGOs (ILIGO), Enhanced LIGOs (ELIGO), ALIGOs 
and ETs. 


Table 2 The S/N for pairs of ILIGO, ELIGO, ALIGO and 
ET 


System 

ILIGO 

(S/N) 

ELIGO 

ALIGO 

ET 

DNS 

3.1 X 10-2 

2.7 

12.0 

1.6 X 10^ 

BHNS 

5.6 X 10-2 

4.2 

16.0 

3.1 X 10^ 

BBH 

2.0 X 10-3 

0.18 

0.69 

1.9 X 103 


For a cross correlation of two interferometric de¬ 
tectors, we can quantify the detectability of a stochas¬ 
tic background by calculating the signal-to-noise ratio 
(S/N), which in this case is given bv[55]: 


(S/Nf 





(54) 


where 5/ and S'^ are the spectral noise densities, T 
is the integration time (we are considering T = lyr), 
7 (z^) is the overlap reduction function, which depends 
on the relative positions and orientations of the two in¬ 
terferometers; and I7gw is the energy density parameter, 
which is given bv[4()) 

- (55) 






BG 


In Table [2] we show the results for S/N for the three 
families of compact binaries. For pairs of ILIGO, we 
have S/N< 1 for the three families of compact binaries, 
which indicates that there is no possibility of detec¬ 
tion; for pairs of ELIGO and ALIGO we have S/N> 1 
for the BHNS and DNS systems, where BHNS systems 
have the highest possibility of detection; for pairs of 
ET, the S/N is higher than ^ 10^ for the three families 
of systems. 


7 Conclusions 

Our main aim here was to develop an alternative method 
to calculate the background generated by cosmologi¬ 
cal compact binary systems during their periodic or 
quasiperiodic phases. We use here an analogy with a 
problem of statistical mechanics in order to perform 
such a calculation as well as taking into account the 
temporal variation of the orbital parameters of the sys¬ 
tems. Such a background is of particular interest since it 
could well form an important foreground for the planned 
gravitational wave interferometers DEGIGO, BBO, LISA 
or eLISA, ALIGO and ET. 

This new tool for the calculation of the background 
has the advantage of being simple and versatile, because 
the distribution functions and other functions and pa¬ 
rameters could be easily changed without the need of 
modifying the formalism. 
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Our results show that the background generated by 
the cosmological compact binaries during their periodic 
and quasiperiodic phases does not form a foreground 
for LISA, eLISA, BBO, DECIGO and ALIGO. For the 
ET a foreground (confusion noise) could marginally be 
formed. It is worth recalling that in the present study, 
we only considered systems in circular orbits, and the 
inclusion of the eccentricity will be investigated in a fu¬ 
ture study. Besides, we will use the formalism developed 
here in the calculation of the background generated by 
the coalescence compact binary systems. 

We also investigated if the backgrounds studied here 
could be detected using a cross-correlation of two inter¬ 
ferometric detectors. We note that pairs of ALIGOs and 
ETs could in principle detect such backgrounds. 

Given that by the time ET becomes operative, there 
will probably be some advanced configuration of VIRGO, 
and then, it will be possible to correlate them. Such 
a cross correlation would be, in principle, much more 
sensitive than a pair of ALIGOs. Therefore, it will be 
possible to investigate the kind of backgrounds studied 
here. 
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A Appendix 


Using the definition given by 1(23, we have 

dn = H{u)du, (56) 


which is the fraction of systems originated at time to with 
orbital frequencies in the interval dv. Using the binary for¬ 
mation rate corresponding to each type of system, we have 

(57) 

where the subindex bin refers to each type of system (DNS, 
BHNS and BBH). Now, integrating over dto, we get 


dn 

dV 


j 'ri\,i^{to)H{v)dto 


dv, 


(58) 


where the expression in brackets is the number of systems per 
unit frequency interval and per comoving volume, which may 
be used as a distribution function in the calculation of the 
background. 
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